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Preface 

The  purpose  of  this  study  was  to  develop  a  method  for 
predicting  realistically  curved  fallout  contours  with  variable 
winds  in  smearing  codes.  The  method  requires  two  steps: 

(1)  location  of  the  curved  hotline,  and  (2)  calculation  of 
off-axis  dose  rates.  A  curved  hotline  locator  model  was 
developed  from  correlated  DELFIC  wafer  height  data  and  a 
real  wind  particle  translation  scheme.  The  model  was  tested 
by  comparing  results  to  DELFIC  output  and  to  Castle- Bravo 
fallout  footprints.  The  two-dimensional  dose  rate  integral 
was  solved  analytically,  enabling  the  calculation  of  off- 
axis  dose  rates. 

This  two-step  method  can  be  adapted  to  existing  smearing 
codes  by  incorporating  the  hotline  locator  model  and  the 
dose-rate  equation. 

I  thank  Dr.  C.J.  Bridgman  for  his  patience  and  guidance 
during  this  thesis  study.  In  addition,  I  thank  Major  W.S. 
Bigelow  for  providing  the  microfiche  DELFIC  data  that  I  used 
for  the  stabilized  cloud  fits,  and  the  Defense  Nuclear  Agency 
for  providing  Castle-Bravo  documentation  and  movies.  Finally, 
thank  you,  Jeannie. 

A.T.  Hopkins 
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'  Abstract 

A  method  was  developed  to  treat  non-constant  winds  in 
fallout  smearing  codes.  The  method  consists  of  two  steps: 

(1)  location  of  the  curved  hotline 

(2)  determination  of  the  off-hotline  activity. 

To  locate  the  curved  hotline,  the  method  begins  with  an 
initial  cloud  of  20  discretely-sized  pancake  clouds, 
located  at  altitudes  determined  by  weapon  yield.  Next, 
the  particles  are  tracked  through  a  300  layer  atmosphere, 
translating  with  different  winds  in  each  layer.  The  connec¬ 
tion  of  the  20  particles'  impact  points  is  the  fallout 
hotline.  The  hotline  location  was  found  to  be  independent 
of  the  assumed  particle  size  distribution  in  the  stabilized 
cloud.  The  off-hotline  activity  distribution  is  represented 
as  a  two-dimensional  gaussian  function,  centered  on  the 
curved  hotline.  Hotline  locator  model  results  were  compared 
to  numerical  calculations  of  a  hypothetical  100  kt  burst 
and  to  the  actual  hotline  produced  by  the  Castle-Bravo  15  Mt 
nuclear  test. 


A  TWO  STEP  METHOD 


TO  TREAT  VARIABLE  WINDS 
IN  FALLOUT  SMEARING  CODES 

I .  Introduction 

Background 

An  atmospheric  nuclear  burst  produces  a  cloud  of 
radioactive  particles  which  translate  with  the  winds  as 
they  fall.  Military  and  civilian  studies  (Refs  1  through 
4)  use  both  numerical  and  analytic  fallout  prediction 
models  to  compute  fallout  activity  levels  as  a  function  of 
location.  Analytic  models  (Refs  1  and  2)  are  used  for  most 
studies,  primarily  because  they  are  economical.  Economy 
is  achieved,  in  part,  by  using  a  constant  wind  vector 
(speed  and  direction)  to  smear  radioactivity  along  the 
ground.  However,  the  constant  wind  assumption  is  a  severe 
restriction  on  the  smearing  model  which  is  universally 
recognized  (Ref  1).  The  smearing  code's  ground  trace  is 
symmetric  about  the  downwind  axis,  or  hotline,  the  locus  of 
peak  activity  levels.  These  hotlines  are  always  straight 
because  the  wind  is  assumed  constant  at  all  altitudes 
beneath  the  falling  cloud.  Figure  1  shows  fallout  contours 
obtained  from  constant  vs  variable  winds  in  a  smearing 
calculation . 
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Figure  1.  Fallout  Contours 


Problem  Statement 


Fallout  smearing  codes  need  to  account  for  spatially 
and  temporally  varying  winds.  Specifically,  smearing  codes 
need  techniques  to  locate  curved  hotlines  and  to  compute 
off-axis  dose  rates. 

Approach 

Two  steps  are  taken  to  model  variable  winds  in  fallout 
smearing  codes.  First,  the  hotline  is  located.  Second,  the 
off-hotline  dose  is  predicted. 

To  locate  the  hotline,  a  model  was  developed  to  compute 
trajectories  of  20  different  particle  sizes.  The  particles 
start  from  altitudes  computed  with  empirical  fits  to  numer¬ 
ical  code  data,  and  fall  through  a  discretely-layered 
atmosphere.  Each  atmosphere  layer  is  100  meters  thick  and 
the  model  has  300  layers.  The  Davies -McDonald  method  was 
used  to  compute  terminal  velocities  in  model  layers.  The 
horizontal  motion  of  each  particle  is  computed  and  summed 
from  starting  altitude  to  the  ground.  The  fallout  hotline 
is  the  connection  of  the  20  particles'  landing  points. 

To  predict  dose  rates,  the  two-dimensional  dose  rate 
integral  was  solved  analytically.  The  solution  will  enable 
the  smearing  codes  to  calculate  dose  rates  as  a  function  of 
two-dimensional  winds. 
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Organization 

Section  II  briefly  describes  the  two  principal  types 
of  fallout  codes.  Sections  III  through  VII  discuss  the 
first  step:  hotline  location.  Section  III  quantifies  the 
working  hypothesis  of  a  stabilized  cloud  with  particle 
sizes  that  have  been  gravity-sorted  in  altitude.  Section  IV 
describes  the  hotline  locator  model  with  altitude-varying 
winds.  Sections  V  and  VI  contain  results  of  case  studies  in 
which  computed  hotlines  are  compared  to  numerical  model 
results  and  an  actual  nuclear  test  hotline.  Section  VII 
describes  a  method  to  account  for  temporal  variations  in 
winds . 

Section  VIII  presents  the  second  step:  the  analytic 
solution  of  the  two-dimensional  dose  rate  integral. 

The  Appendices  contain  detailed  descriptions  of  the 


atmosphere  equations  and  the  hotline  locator  model. 


II .  Current  Fallout  Prediction  Models 

This  section  briefly  describes  the  two  principal 
types  of  fallout  codes:  disc  tosser  and  smearing. 

Disc  Tosser  (DELFIC) 

The  Department  of  Defense  Land  Fallout  Interpretive 
Code  (DELFIC)  is  the  standard  numerical  fallout  code  for 
military  studies  (Ref  1).  DELFIC  is  a  multi-dimensional 
model  which  describes  the  motions  of  fallout  particles 
contained  in  discrete  pancake  wafers  of  discrete  particle 
size  groups.  The  DELFIC  data  (Rei  5)  used  in  this  study 
represented  the  particle  clouds  with  one  hundred  discrete 
particle  size  groups  (3  pm  to  4.6  mm)  and  at  least  twenty 
wafers  per  size  group.  Wafer  motion  is  modeled  in  DELFIC 
as  the  fireball  rises,  cools  and  expands  and  as  the  stabi¬ 
lized  cloud  falls  through  real  winds.  As  a  research  tool, 
DELFIC  provided  valuable  data  for  empirical  fits  of  stabi¬ 
lized  wafer  heights  and  for  benchmark  comparisons  of  curved 
hotlines.  However,  DELFIC  is  too  costly  and  detailed  to 
use  when  large  numbers  of  hypothetical  fallout  cases  must 
be  analyzed.  Fast-running,  economical  analytic  codes 
(Refs  1  and  2)  have  been  developed  to  support  parametric 
studies , 


Smearing  Codes  (AFIT  and  WSEG 


The  AFIT  (Air  Force  Institute  of  Technology)  and 
WSEG  (Army  Weapon  Systems  Evaluation  Group)  codes  (Refs  1 
and  2)  are  examples  of  fallout  smearing  prediction  models. 

Both  AFIT  and  WSEG  models  simulate  particle  transport  and 
deposition  by  smearing  the  falling  radioactive  cloud  along 
the  ground,  as  a  constant  wind  vector  translates  airborne 
particles  downwind.  The  radioactive  fallout  is  symmetrically 
distributed  on  the  ground,  normal  to  the  downwind  axis  (hot¬ 
line)  of  the  fallout  footprint.  Smearing  codes  can  only 
predict  unidirectional  hotlines  because  they  assume  a  constant 
wind  vector  for  particle  transport.  The  fundamental  equation 
for  computing  ground  dose  rate  by  smearing  the  descending 
cloud  is:  (Ref  1) 


6(x,y)  =  KYk  /  f  (x,y,t')g(t')dt' 
0 


(II-l) 


D  =  Dose  rate  at  ground  coordinates  (x,y)  R°hbur 


K  =  Source  normalization  constant 

.  ,,  n-rm  Roentgens  miles2  A 

(typically,  2350  Hlotohs  } 


Yield  (kilotons,  kt) 


f(x,y,t')  *  normalized  cloud  spatial  distribution 

(miles’ 2 ) 
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HiiTriT  i.Vii.  i . ‘  '  - 


g(t) 


time-dependent  arrival  rate  of 
radioactivity  on  the  ground  (hour-1) 


x  =  distance  downwind  (mi) 
y  =  distance  crosswind  (mi) 


The  normalized  spatial  distribution  function  describes  the 
lateral  (crosswind  and  downwind)  extent  of  the  radioactive 
ground  trace.  Assuming  gaussian  distributions  for  the 
x  and  y  dependence  of  f(x,y,t')  ,  the  activity  function  is: 


f(x,y,t')  = 


/n Fax  (t^) 


EXP 


-h(- 


x-v  t' 


/TirOy  (t') 


EXP 


'  y-v  t'~ 
-ht—y-) 
°Y  J 


(II- 


=  standard  deviation  in  x-direction  (mi) 

Oy  =  standard  deviation  in  y-direction  (mi) 

vx  =  wind  velocity  in  x-direction  (mi/hr) 

Vy  =  wind  velocity  in  y-direction  (mi/hr) 

Both  smearing  codes  solve  the  dose  rate  integral  by  assuming 
an  "effective"  constant  wind  (Ref  8)  which  is  actually  a 
particle  residence  time-weighted  average  of  a  vertical  wind 
profile.  Essentially,  the  x-direction  is  defined  as  downwind 


% 


and  v  is  zero.  The  dose  rate  integral  then  contains 
a  single  time -dependent  exponential  term  and  is  solved 
analytically: 


D 


SH  vxoy  Cta) 


°Y 


(II ~3) 


Equation  ( 1 1 - 3 )  is  the  working  equation  for  dose  rate 
in  both  the  WSEG  and  AFIT  smearing  codes  (Refs  1  and  2). 
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III.  The  Initial  Stabilized  Cloud 

Description 

A  surface  nuclear  burst  produces  a  fireball  which  rises, 
expands  and  cools,  entraining  air  and  surface  debris.  The 
buoyant  cloud  contains  vaporized,  radioactive  particles  which 
eventually  recondense  on  microscopic  nuclei  or  on  the  sur¬ 
faces  of  other  particles.  When  the  gravitational  force  on 
a  particle  exceeds  the  force  induced  by  updraft  in  the  rising 
cloud,  the  particle  falls  from  the  cloud.  Therefore,  larger, 
heavier  particles  should  begin  to  fall  from  the  rising  cloud 
sooner  than  small  particles.  However,  sorting  by  terminal 
velocity  is  not  perfect  because  particles  have  different 
shapes  and  densities,  and  because  the  flow  inside  the  rising 
cloud  is  not  a  simple  updraft.  Toroidal  flow  and  temperature- 
dependent  changes  in  the  air  properties  help  to  mix  the  par¬ 
ticles  during  cloud  rise  to  stabilization. 

Pancake  Cloud 

The  smearing  codes  (Refs  1  and  2)  assume  that  all  fallout 
particles  fall  from  altitudes  related  to  the  stabilized  cloud 
center  height.  The  WSEG  model  assumes  that  radioactivity  is 
normally  distributed  about  a  fallout  cloud  altitude  given 
by:  (Ref  2) 

HC  *  44.0+6.1  •  AnY- .205 (ZnY+2 . 42) | AnY+2. 42 |  (III-l) 
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HC 


Cloud  Height  (kilofeet) 


Y  *  Weapon  Yield  (megatons,  Mt) 

Bridgman  (Ref  1)  showed  that  the  vertical  cloud  distribution 
can  be  reduced  to  a  pancake  cloud  at  HC.  In  fact,  Reference 
1  demonstrates  that  the  pancake  cloud  is  a  convenient,  accurate 
simplification  of  the  nuclear  fallout  cloud  if  the  vertical 
distribution  of  activity  is  symmetric  about  HC.  However,  when 
particles  start  from  pancake  clouds  and  fall  through  real, 
variable  winds,  they  impact  in  a  nearly  straight  line.  There¬ 
fore,  the  AFIT  simplification  of  a  pancake  initial  cloud 
cannot  be  used  to  predict  curved  hotlines  with  real  winds. 

Distributed  Cloud 

DELFIC  data  (Ref  5)  showed  that  the  particles  in  stabi¬ 
lized  clouds  are  gravity  sorted,  with  large  particles 
generally  lower  in  altitude  than  small  particles.  To  quan¬ 
tify  the  distributed  initial  cloud,  the  DELFIC  data  (Ref  5) 
was  used  to  generate  a  generalized  initial  condition. 

DELFIC' s  Advective  Transport  Module  divides  an  assumed 
particle  size  distribution  into  100  discrete  groups  and 
distributes  each  group  vertically  into  approximately  20  wafers. 
Each  wafer  has  finite  dimensions.  A  total  of  2348  wafer  top 
and  bottom  heights  were  used  for  this  study.  For  each  yield 
(15  Mt ,  1  Mt ,  100  kt ,  10  kt ,  1  kt)  and  particle  size,  wafer 


height  was  taken  to  be  the  middle  of  each  wafer  top  and  wafer 
bottom  height.  Figure  2  shows  a  typical  range  of  wafer 
heights  for  particles  in  a  stabilized  cloud. 

Linear  Regressions  of  Wafer  Center  Heights 

To  obtain  a  single-valued  function  of  particle  diameter 
vs  height  for  use  in  the  smearing  codes,  the  spatial  center 
height  for  each  size  wafer  set  was  computed  as  the  average 
of  the  wafer  center  heights  for  each  particle  diameter. 

Plots  of  the  spatial  center  heights  vs  particle  diameters 
for  each  yield  showed  that  a  linear  fit  accurately  represented 
the  vertical  distribution.  If  the  average  wafer  height  of 
each  size  represents  the  maximum  radioactivity  for  that  size, 
then  the  locus  of  points  along  the  ground,  where  the  peak  of 
each  size  falls,  would  be  the  fallout  hotline.  To  test  the 
"peak  radioactivity  height"  hypothesis,  this  thesis  study 
includes  a  direct  comparison  with  the  DELFIC  hotline  for  a 
100  kt  surface  burst.  The  hotlines  match  very  well,  suggest¬ 
ing  that  size-dependent  radioactivity  peaks  can  be  modeled 
at  spatial  average  wafer  heights  in  the  stabilized  cloud. 

The  case  study  is  described  in  Section  V.  The  five  DELFIC 
output  listings  (1  kt ,  10  kt ,  100  kt ,  1  Mt ,  15  Mt)  provided 
data  to  derive  five  linear  fits.  Particle  diameters  were 
correlated  with  spatial  central  wafer  heights;  the  data  and 
lines  are  shown  in  Figure  3.  Higher  yields  generate  larger 
vertical  "spreads,"  steeper  curves  and  higher  clouds. 
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Figure  2.  Typical  Range  of  DELFIC  Wafer  Heights 


in  a  Stabilized  Cloud 


Following  are  the  equations  developed  from  linear  regres¬ 
sions  : 


H  =  (SLOPE) (D  )  +  INTERCEPT 

c  p 


Hc  =  Average  height  of  wafer  centers 

in  stabilized  cloud  (m) 


Dp  .  =  Particle  diameter  (micrometers) 


SLOPE 

=  (meters/micrometer) 

INTERCEPT 

=  Altitude  of  D  = 
P 

Linear  Fit  (m) 

0  on  a 

Y 

SLOPE 

INTERCEPT 

1 

kt 

-4.82 

2670  (m) 

10 

kt 

-5.44 

5581 

100 

kt 

-7.22 

9496 

1 

Mt 

-10.  2 

13677 

15 

Mt 

-17.4 

26643 

Polynomial  Least  Squares  Fits 

To  obtain  similar  relations  for  other  yields,  the 
slopes  and  intercepts  were  correlated,  using  a  polynomial 
least -squares  program  (Ref  12).  The  following  fits  were 
obtained : 


(III. 2) 
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SLOPE  «  -EXP(1.S74-. 01197 i.nY+ .03636 (AnY) 2 

-  .  0041  ( ilnY)  3+.000l965(  £nY)  4  )  (III-3) 

INTERCEPT  =  EXP(7.889+.  3  4S,nY+  .001226  (£nY)  2 

-  .  005227  (£nY)  3+.000417  (£nY)‘*)  (II 1-4) 

Y  =  weapon  yield  (kilotons,  kt) 

Figure  4  shows  the  variations  in  slopes  (absolute  values) 
and  intercepts  with  weapon  yield.  Given  a  weapon  yield, 
the  polynomial  fits  are  used  to  generate  values  of  slopes 
and  intercepts,  defining  the  linear  equation  for  vertical 
distribution  in  a  stabilized  cloud.  So,  in  a  stabilized 
cloud,  the  central  wafer  height  of  a  particle  with  diameter, 
d(  m)  ,  is  given  by  (III-2).  Figure  5  shows  the  parametric 
variation  of  central  wafer  height  with  weapon  yield  for 
thirteen  particle  sizes.  As  expected,  larger  yields  raise 
wafers  higher  in  the  atmosphere. 

Summary 

The  assumption  of  a  single  pancake  initial  cloud  used 
in  the  AFIT  model  has  been  replaced  with  a  vertical  array 
of  pancake  clouds,  one  for  each  particle  size.  Particle  size 
vs  height  in  the  initial  cloud  is  computed  from  an  empirical 
fit  to  DELFIC  data. 
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IV.  Hotline  Locator  Model 


Fallout  hotlines  curve  realistically  in  DELFIC  calcu¬ 
lations  because  DELFIC  traces  the  motions  of  wafers  contain¬ 
ing  single-size  particles  through  real  wind  profiles.  The 
smearing  codes  compute  uni-directional  fallout  contours  with 
a  constant  "effective"  wind  described  in  References  2  and 
8.  However,  real  winds  can  vary  both  (spatially)  throughout 
the  altitudes  where  particles  fall  and  (temporally)  during 
the  time  that  the  cloud  resides  in  the  atmosphere.  This 
section  describes  the  computer  model  which  predicts  hotline 
locations  with  spatial  wind  variations  beneath  the  stabilized 
cloud.  Section  VII  describes  a  technique  to  simulate  the 
combined  effects  of  spatial  and  temporal  variations. 

Description 

The  hotline  locator  model  computes  and  summarizes  the 
motion  of  the  initial  cloud  described  in  Section  III.  The 
model  contains  an  atmosphere  which  is  divided  into  three- 
hundred  100-meter-thick  layers,  with  state  properties  com¬ 
puted  as  described  in  Appendix  A.  Figure  6  shows  the  model 
structure.  With  a  real  wind  profile  as  input,  the  wind 
velocity  and  direction  are  linearly  interpolated  between 
data  points,  and  a  two-component  wind  vector  (East:  positive 
x,  North:  positive  y)  is  generated  for  each  atmospheric 
layer.  Twenty  particles,  representing  a  fallout  size 
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spectrum,  fall  from  the  initial  altitudes.  Altitudes  are 
computed  for  a  gravity-sorted  stabilized  cloud  using  poly¬ 
nomial  fits  and  linear  height  functions.  Particle  sizes 
were  selected  to  represent  the  fallout  population  and  to  give 
a  downrange  spread  (hotline  length)  comparable  with  published 
data  (Refs  5  and  14).  Twenty  particle  diameters  are  speci¬ 
fied  and  their  impact  points  (x,y)  and  arrival  times  are 
determined  by  winds,  atmospheric  state  variables  and  particle 
residence  times  in  each  model  layer. 

Davies -McDonald  Method 

Particle-layer  residence  times  were  computed  from 
terminal  velocities  which  were  calculated  with  the  Davies 
McDonald  method  (Refs  6  and  7).  Davies  showed  that  Reynolds 
number.  Re  ,  of  a  moving  sphere  is  related  to  the  product 
of  drag  coefficient,  ,  and  Re2  ,  and  the  functional 

dependences  were  published  in  Reference  7: 

If  CpRe2  <  140  or  Re  <  4 

Re  =  CDRe2  -  2 . 3363x10'“  (CpRe2 ) 2 


+  2. 0154x10' 6  (CjjRe2 ) 3  -  6 . 9104x10' 9 (CpRe2) 4  (IV-1) 


If  100  <  CDRe2  <  4 . 5x10 7  or  3  <  Re  <  1000 
LoglQRe  =  -1.29536  +  .  986  (log^CpRe2 ) 

-  .046677*(logl0CDRe2) 2 

+  .0011235(log10CDRe2)3 

McDonald  (Ref  6)  showed  that  CQRe2  can  be  obtained 
explicitly  for  a  given  particle  size  and  density  at  an 
altitude  with  known  kinematic  viscosity  and  density: 

CDRe2  =  8W/irpv2  , 


where 


Cn  =  drag  coefficient 


'D 

Re 


=  Reynolds  number  = 


Vzd 


p  =  density  of  air  (g/cm3) 

Vz  =  terminal  velocity  of  particle  (cm/sec) 

d  =  particle  diameter  (cm) 

v  *  kinematic  viscosity  of  air  (cm2/sec) 

W  *  weight  of  particle  ■  yirCj)  3p  g  (-%—-) 

sec 2 
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Ps  =  density  of  solid  particle  (g/sec2) 

g  =  gravitational  constant  (cm/sec2) 

Centimeter-gram-second  (cgs)  units  were  used  to  facilitate 
comparisons  with  data  published  in  References  6  and  7. 
Terminal  velocity  of  a  particle  at  a  specified  altitude  is 
obtained  by  (a)  computing  C^Re2,  (b)  using  Davies'  equations 
to  obtain  Re,  and  (c)  from  the  definition  of  Re: 


V 


z 


Rev 

T" 


The  residence  time,  x  (sec),  of  a  particle  in  a  model  layer 
of  thickness  AZ  was  computed  from  the  average  velocity  of 
the  particle  in  the  layer. 


AZ 


Tik 

L  2 

where 

Tik 

residence  time  of  it^1  particle  in 
layer  k  (sec) 

AZ  * 

thickness  of  layer  (cm) 

- 

terminal  velocity  of  particle  at  bottom 
height  of  layer,  level  j  (cm/sec) 

<Vj.i  a 

terminal  velocity  of  particle  at  top 
height  of  layer,  level  j+1  (cm/sec) 

(layer  k  is  bounded  by  heights  j  and  j+1) 
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Knowing  the  residence  time,  particle  translation  in  each 
layer  was  computed  from  the  layer's  mean  wind  components. 
Hotline  coordinates  (x,y,  on  the  ground)  were  obtained  by 
summing  translations  in  all  layers. 


£  Vxk  Tik 
HP 


l  Vyk  Tik 


where 


x  =  ground  coordinate  in  x  (east)  direction  (m) 

y  =  ground  coordinate  in  y  (north)  direction  (m) 

{jc  3  x  component  of  wind  velocity  in  layer  k  (m/sec) 

^  -  y  component  of  wind  velocity  in  layer  k  (m/sec) 

■Ly  =  residence  time  of  it^1  particle  size  in 
layer  k  (sec) 


H  =  height  of  ground,  MSL  (m) 
§ 


Hp  =  height  of  particle  in  stabilized  cloud  (m) 


Summary 


The  model  computes  the  total  horizontal  movement  of 
20  particles  as  they  fall  from  an  initial,  vertically- 
distributed  cloud.  Real  wind  data  is  interpolated  into 
each  model  layer.  Particle  fall  speeds  are  determined  by 


L 
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the  Davies -McDonald  method.  The  hotline  is  the  connection 
of  the  particles'  ground  impact  points.  Appendix  B  des¬ 
cribes  the  computational  algorithm  and  includes  a  code 
i  listing. 


V.  Case  b. ndy  1:  Washington  DC, 

3  Oct  68,  100  kt  Hypothetical 

In  Reference  4,  Norment  compared  the  fallout  contours 
generated  by  DELFIC,  WSEG  and  SEER  in  real  winds  using  a 
set  of  hypothetical  and  real  atmospheric  nuclear  detonations. 
To  compare  hotline  locations  generated  in  this  thesis  study 
with  DELFIC  results,  the  Washington  DC,  3  Oct  68,  100  kt 
case  was  selected,  primarily  because  the  DELFIC  hotline 
exhibited  curvature.  The  wind  profile  is  shown  in  Figure  7. 

DELFIC  Hotline  and  Size  Distribution  Independence 

Dose  rate  contours  from  DELFIC  (Ref  4)  are  shown  in 
Figure  8.  The  double  hotline  was  confirmed  by  analyzing 
the  DELFIC  numerical  output  for  this  case  provided  in 
Reference  5.  Two  peaks  in  dose  rate  were  observed  at  down- 
range  distances  greater  than  -v35  km.  Figure  9  is  a  plot 
of  two  sets  of  DELFIC  hotline  coordinates,  generated  by 
two  different  particle  size  distributions  (PSD).  As  shown 
in  Table  I,  the  default  PSD  contains  an  average  particle 
diameter  of  .407  pm,  and  the  large  PSD  contains  a  24.0  pm 
average  diameter.  In  the  region  of  the  double  hotline,  the 
bottom  line  has  higher  dose  rates,  suggesting  that  the 
hotline  is  curving  downward  as  the  cloud  translates  down¬ 
wind.  The  two  DELFIC  calculations  were  analyzed  to  identify 
downrange  locations  of  peak  dose  rate,  thereby  defining 
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Figure  7. 


Wind  Profile,  Washington,  D.C.,  3  Oct  68 
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Figure  8.  DELFIC  Dose  Rate  Contours  for  Wash.,  D.C., 
3  Oct  68,  lOOkt  Case  (Norment) 
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Figure  9.  Hotlines  (DELFIC)  Washington,  D.C.,  3  Oct  G8,  lOOkt 
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the  hotline.  Figure  9  also  shows  the  hotline  provided  by 
the  larger  distribution.  The  DELFIC  hotline  locations  are 
identical,  independent  of  the  two  different  particle  size 
distributions  assumed  for  the  calculations.  For  a  given 
yield,  DELFIC  raises  wafers  of  fixed  particle  size  to  the 
same  heights  in  the  stabilized  cloud.  Since  the  two  size 
distributions  overlap,  DELFIC  predicted  the  same  hotlines. 

In  the  region  where  the  distributions  do  not  overlap,  the 
hotlines  are  continuous.  So,  given  a  yield  and  a  wind 
profile,  DELFIC  will  predict  constant  hotline  locations. 
Therefore,  a  hotline  locator  model  should  predict  the  same 
hotline  locations  with  different  assumed  particle  size  dis¬ 
tributions.  In  fact,  the  hotline  locator  model  contains 
an  array  of  particle  sizes,  selected  only  to  resolve  hotline 
coordinates . 

Hotline  Comparisons 

The  hotline  locator  model,  described  in  Section  IV, 
was  used  to  predict  the  Washington  DC  100  kt  hotline. 
Particle  sizes  selected  for  the  calculations  are  shown  in 
Table  I.  Particle  sizes  were  selected  from  the  DELFIC 
calculation  to  represent  (approximately)  intervals  of  5% 
in  cloud  activity.  The  wind  profile  from  Figure  7  was  used 
as  input. 
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Table  I 

Particle  Size  Distributions 


N  (d)  =  - i -  EXP-Ss 

/!¥  $,n(r)d  L  n° 


N(d)  =  Number  per  micrometer  size  range 
d  =  Particle  diameter  (m) 


Default  PSD 
Large  PSD 


Median  Diameter 
(T,  ym 
0.407 


Deviation 
o ,  urn 
24.0 


Table  II 

Particle  Size  Arra 


.  0605 

.  0089 

.  0290 

.  0081 

.0218 

.  0074 

.0184 

.0067 

.0160 

.0050 

.0144 

.  0054 

.0129 

.0048 

.0117 

.0042 

.  0106 

.0035 

.  0097 

.  0026 
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Figure  10  shows  two  predictions  from  the  hotline 
locator  model  and  the  DELFIC  hotline.  The  pancake  cloud 
hotline  was  generated  by  starting  all  particles  at  9.2  km 
altitude,  computed  with  Eq  (III-l).  The  wafer  centers 
hotline  (labeled  "wafer  ctrs")  was  produced  by  starting 
the  particles  from  altitudes  computed  from  linear  size- 
height  relations  described  in  Section  III.  All  computa¬ 
tions  used  real  wind  profiles  beneath  the  stabilized  cloud. 

The  pancake  cloud  hotline  is  nearly  unidirectional. 
Computations  performed  for  this  thesis  study  indicate 
that  the  pancake  cloud  assumption  severely  restricts  hot¬ 
line  curvature.  If  all  particles  fall  from  the  same 
altitude,  hotline  curvature  must  be  caused  by  different 
fall  velocities  (residence  times)  in  wind  layers.  From  a 
pancake  cloud,  the  differences  are  not  large  enough  to 
produce  significant  curvature  in  the  particles'  ground 
traces.  The  pancake  cloud  hotline  deviated  (crossrange) 
from  the  DELFIC  hotline  by  approximately  10  km  at  50  km 
downrange.  At  long  distance  (>  200  km)  the  pancake  and 
"wafer  ctrs"  hotlines  agree  within  5  km. 

The  "wafer  ctrs"  hotline  is  on  the  DELFIC  hotline  at 
ranges  less  than  40  km,  where  it  properly  curves  toward 
the  more  intense  part  of  DELFIC' s  dual  ground  trace.  The 
agreement  between  the  "wafer  ctrs"  hotline  and  DELFIC  is 
expected  because  the  "wafer  ctrs"  inputs  are  DELFIC  wafer 
heights.  However,  since  a  single  height  was  used  for  each 
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particle  size  in  the  hotline  locator  model,  this  case 
study  showed  that  the  hotline  location  can  be  computed 
from  a  stabilized  cloud  modeled  with  a  linear  relationship 
between  particle  size  and  height.  Furthermore,  since  the 
hotline  represents  the  locus  of  peak  activity  on  the 
ground,  it  can  be  inferred  that  the  wafer  central  heights 
represent  maxima  in  vertical  activity  distributions  for 
corresponding  particle  sizes  in  the  initial  stabilized 
cloud. 

Sensitivities 

For  100  kt,  Eqs  CHI-3)  and  (III-4)  give  7.228  JjSL 
and  9496  m  for  slope  and  intercept  of  the  particle  size 
vs  height  array  in  stabilized  clouds.  To  evaluate  the 
sensitivity  of  hotline  location  to  consistent  variations 
in  the  vertical  array,  additional  calculations  were  per¬ 
formed  with  yield  varying  by  ±  10%  and  ±  20%.  Sensitivity 
of  stabilized  cloud  height  is  shown  in  Figure  11,  in  terms 
of  variations  in  the  mean  wafer  center  heights.  A  variation 
of  ±  20%  in  weapon  yield  produced  about  ±  5%  variation  in 
height  over  the  entire  range  of  particle  sizes.  Hotline 
location  was  relatively  insensitive  to  variations  (uncer¬ 
tainties)  in  weapon  yield.  Figure  12  shows  the  maximum 
hotline  excursion,  obtained  for  Yield  =  80  kt . 

In  summary,  the  Washington  DC  case  study  supports 
the  following  conclusions: 
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Id  Sensitivity  of  the  Stabilized  Cloud 


Figure  12.  Hotlines,  80kt  Sensitivity 


1.  DELFIC  hotline  location  is  independent  of  the 
assumed  particle  size  distribution; 

2.  The  hotline  location  can  be  computed  from  a 
linear  relation  between  particle  size  and 
height  in  the  initial  stabilized  cloud. 
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VI.  Case  Study  2:  Castle-Bravo 

4  Mar  54,  15  Mt ,  Pacific  Test 

The  hotline  locator  model  was  used  to  compute  the 
hotline  coordinates  of  the  Castle-Bravo  nuclear  cloud. 
Castle-Bravo  was  also  predicted  by  DELFIC,  so  this  case 
study  was  an  opportunity  to  compare  the  hotline  locator 
results  to  both  DELFIC  and  actual  fallout  contours.  When 
it  was  learned  that  the  models  agreed,  but  neither  pre¬ 
dicted  the  actual  hotline,  a  new  technique  was  developed 
to  back  calculate  initial  particle  heights  in  the  cloud. 
These  back  calculations  provided  a  new  interpretation  of 
the  stabilized  cloud  produced  by  the  high  yield  Pacific 
test . 

AFSWP  Fallout  Contours 

The  Castle-Bravo  nuclear  weapon  test  was  detonated 
in  the  Bikini  Atoll  on  4  March  1954.  The  unexpectedly 
high  fission  yield  and  changes  in  the  wind  structure 
caused  significant  fallout  to  reach  populated  islands 
(Ref  13).  The  Armed  Forces  Special  Weapons  Project  (AFSWP) 
collected  air  and  ground  fallout  samples  and  constructed 
the  fallout  contour  map  shown  in  Figure  13  (Ref  14).  The 
fallout  hotline  dipped  southward  before  tracking  to  the 
east  and  northeast.  Figure  13  shows  estimates  of  fallout 
arrival  times  at  various  distances  downrange.  Wind  data 


37 


was  obtained  from  Enewetok,  Kwajalein  and  Rongerik  Islands, 
as  well  as  the  U.S.  Navy  ship  S.S.  Curtis.  The  best  fall¬ 
out  predictions  were  calculated  with  a  wind  profile  devel¬ 
oped  from  data  reported  by  Dean  and  Olmstead  (Ref  4),  and 
plotted  in  Figure  14. 

Back  Calculations 

In  this  study,  the  hotline  locator  model  was  used  to 
predict  the  Bravo  hotline.  The  prediction  generally 
agreed  with  DELFIC,  but  did  not  agree  with  the  experimental 
(AFSWP)  hotline.  Figure  15  shows  the  DELFIC  prediction 
(Ref  4).  Figure  16  shows  the  AFSWP  hotline  and  two  pre¬ 
dictions  with  the  hotline  locator  model.  "Wafer  ctrs" 
refer  to  calculations  with  the  particle  starting  heights 
computed  from  polynomial  correlations  of  linear  fit  para¬ 
meters.  The  pancake  hotline  was  produced  by  a  single  pan¬ 
cake  cloud  located  at  16.9  km  (computed  with  Eq  (III-l). 

Dean  and  Olmstead  winds  were  used  to  translate  the  particles. 

To  help  explain  the  Bravo  fallout  motion,  this  study 
used  real  winds  in  a  "vertical  stack"  diagnostic  technique 
to  back  calculate  the  initial  heights  in  the  stabilized 
cloud  from  which  particles  had  to  fall  in  order  to  have 
landed  on  the  AFSWP  hotline.  The  "vertical  stack"  technique 
computes  ground  traces  of  continuous  single-size-particle 
stacks  from  surface  to  30  km  altitude.  The  intersections 
of  vertical  stack  ground  traces  with  the  actual  hotline 
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Figure  15.  DELFIC  Calculation  of  Castle-Bravo 
Fallout  Contours  (Norment) 
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Figure  16.  Hotlines,  Castle-Bravo  (Pacific)  1  Mar  54  15Mt 

Dean  and  Olmstead  T-0  Winds 


define  the  starting  altitudes  of  specific  particle  sizes. 
Different  size  particles  have  different  ground  traces 
because  of  differences  in  terminal  velocities,  residence 
times,  and  vector  translations  in  wind  layers. 

Figures  17  and  18  show  ground  traces  of  vertical 
stacks.  In  Figure  17,  the  stacks  extended  up  to  the 
lowest  wafer  bottom  height  given  for  each  particle  size 
in  the  corresponding  DELFIC  calculation  (Ref  5).  Lowest 
wafer  bottom  heights  are  shown  in  Table  III.  The  lowest 
DELFIC  wafers  landed  closest  to  the  AFSWP  hotline. 

Figure  17  shows  that  DELFIC  could  not  possibly  have  pre¬ 
dicted  the  AFSWP  hotline,  since  DELFIC 's  stabilized 
cloud  was  too  high  for  wafers  to  land  on  the  hotline. 

In  Figure  18,  the  stacks  extended  up  to  30  km,  the  maximum 
level  height  in  the  hotline  locator  model.  The  back 
calculations  clearly  display  two  features  which  may  explain 
why  previous  predictions  failed  to  replicate  the  south¬ 
ward  curvature  in  the  Castle-Bravo  hotline: 

1.  Within  130  km  of  ground  zero,  particles  which 
landed  on  the  hotline  started  from  altitudes  well 
below  stabilized  cloud  heights  predicted  by 
DELFIC  and  the  WSEG  empirical  equation. 

2.  Between  130  and  280  km,  the  hotline  is  nearly 
parallel  with  particle  stack  ground  traces, 
implying  that  initial  particle  sizes  did  not 
change  much  with  altitude  in  the  portion  of  the 
cloud  that  landed  in  that  region. 
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Figure  17.  Ground  Traces  of  Single-Size  Particle  Stacks, 
up  to  DELFIC  Lovjest  Wafer  Bottom  Heights 


XNiorau 


Figure  18.  Ground  Traces  of  Single-Size  Particle  Stacks,  up  to  30km 


iNiorau 


Figure  19.  Ground  Traces  of  Single-Size  Particle  Stacks, 

up  to  Adjoint  Heights 


Table  III 


DELFIC  IS  Mt  Calculation 
Lowest  Wafer  Bottom  Heights,  Hb 


Diameter  (ym) 

Hb  (km) 

500. 

11.5 

300. 

14.5 

200. 

16.1 

150. 

17.0 

130. 

17.4 

115. 

17.6 

100. 

17.9 

95. 

17.9 

90. 

18.1 

85. 

18.1 

83. 

18.2 

80. 

18.2 

78. 

18.3 

75. 

18.3 

73. 

18.3 

70. 

18.4 

To  define  the  initial  cloud,  stack  heights  were  estimated 
to  produce  ground  traces  that  intersected  the  AFSWP  hotline. 
Figure  19  shows  ground  traces  of  stacks  which  extended  up 
to  heights  given  in  Table  IV.  Each  ground  trace  is  a 
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connection  of  N  points,  where  N  is  the  number  of  layers 
in  a  column  of  fixed  diameter  particles.  The  N  for  each 
particle  size  was  estimated  from  the  fractional  length  of 
the  ground  trace  up  to  the  intersection  with  the  AFSWP 
hotline.  A  first  estimate  for  the  starting  height  of  a 
particle  is: 

L 

H  =  N  x  100  meters  . 
lT 

H  =  height  of  a  particle  stack 
Lu  =  length  of  trace  upwind  of  intersection  point 

LT  =  total  length  of  trace 

N  =  number  of  points  on  trace 

100  meters  =  layer  thickness 

The  method  is  shown  schematically  on  Figure  20.  Fine- 
tuning  is  done  by  trial  and  error,  putting  new  estimates 
of  particle  stack  heights  into  the  model. 

The  actual  Bravo  hotline  (AFSWP)  can  be  approximated 
with  three  segments: 

(1)  x  s  130  km 

(2)  130  km  s  x  s  280  km 

(3)  280  km  s  x 
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North  of  GZ 


Figure  20.  Ground  Trace  Intersection  Method 


Back  calculations  of  the  first  segment  showed  that  particle 
starting  heights  were  inconsistent  with  observed  fallout 
arrival  times.  Furthermore,  the  back  calculated  initial 
cloud  was  not  gravity  size  sorted.  This  observation  sug¬ 
gested  that  the  early  fallout  came  from  a  complex  distri¬ 
bution  of  particle  sizes  at  low  altitudes  (stem).  To 
generate  the  southeastward  trace  of  the  actual  hotline, 
the  first  segment  was  approximated  with  a  pancake  cloud 
located  at  9.2  km  altitude.  The  stack  contours  for  the 
pancake  are  shown  in  Figure  21,  and  the  pancake  hotline 
is  shown  in  Figure  22.  From  Figure  17,  it  was  observed 
that  the  second  segment  is  nearly  identical  to  the  80  pm 
diameter  particle  trace.  Therefore,  the  first  two  segments 
of  the  AFSWP  hotline  (the  curved  part)  can  be  located  by 
computing  fallout  motion  from  a  low  altitude  pancake  cloud 
beneath  a  nearly  single-sized  "mixing  region."  This 
information  and  the  fact  that  some  stack  traces  intersect 
the  hotline  twice,  imply  that  the  Bravo  cloud  had  at  least 
two  distinct  vertical  distributions  of  radioactive  fallout 
particles : 

1.  A  low  altitude  stem  distribution  that  was 
primarily  affected  by  the  wind  profile 
beneath  9.2  km; 

2.  A  well-mixed,  almost  uniform  size  region  that 
accounts  for  the  slow  change  in  particle  size 
along  segment  2. 
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igure  21.  Ground  Traces  of  Single-Size  Particle  Stac 


Figure  22.  Hotlines,  Castle-Bravo,  Pancake 


Figures  23  and  24  show  the  adjoint  hotline  calculation 
with  a  pancake  cloud  to  represent  the  stem  and  a  well- 
mixed  single-size  (80  pm)  layer  between  10  and  16  km 
altitude. 

The  need  for  a  special  treatment  of  the  Castle-Bravo 

fallout  was  recognized  by  Norment  (Ref  4): 

•  .high  yield  explosions  over  coral  and 
ocean  water  need  special  treatment.  .  .Clearly 
a  particle  size  or  size  activity  distribution 
that  is  quite  distinct  from  any  used  here  is 
required  for  this  shot.”  (Ref  4:65) 

Reference  15  further  notes  that  the  WSEG-10  model  produces 

good  agreement  with  the  AFSWP  contours,  "in  spite  of  a 

stabilized  cloud  that  is  half  as  high  as  it  should  be." 

(Ref  15:65)  Table  V  shows  cloud  top  and  bottom  heights 

used  in  Reference  4: 


Table  V 

Castle-Bravo  Cloud  Heights 


Observed 

DELFIC 

WSEG-10 

Cloud 

Base  (km) 

16.9 

20.3 

10.  7 

Cloud 

Top  (km) 

34.7 

39.7 

22.8 

On  the  basis  of  the  back  calculations  for  this  study,  it  can 
be  inferred  that  the  WSEG-10  agreement  is  coincidental, 
caused  by  the  existence  of  a  stem  distribution  which  occurred 
well  below  the  observed  stabilized  cloud  height.  A  different 


Figure  24.  Hotlines,  Castle-Bravo,  9.5km  Pancake  +  Single-Size  Layer 


particle  size  distribution  will  not  improve  code  predictions, 
since  the  DELFIC  hotline  location  is  insensitive  to  size 
distribution.  Rather,  a  different  (stem  +  cloud)  vertical 
spatial  description  of  the  stabilized  cloud  is  needed  to 
explain  the  actual  location  of  the  fallout  hotline. 

Beyond  280  km,  it  is  obvious  from  Figure  18  that  the 
Dean  and  Olmstead  winds  at  burst  time  could  not  have  trans¬ 
lated  particles  onto  the  hotline  segment  3.  In  fact,  using 
shot-time  winds  to  predict  fallout  locations  at  6+  hours 
after  detonation  is  tacitly  assuming  no  temporal  or  spatial 
variations  in  the  wind  profile.  This  assumption  will  intro¬ 
duce  large  uncertainties  in  fallout  location  predictions, 
especially  when  the  distances  traveled  by  the  fallout  parti¬ 
cles  are  on  the  order  of  synoptic  weather  patterns,  or 
10“  -  106  meters  (Ref  15).  Over  a  synoptic  scale,  fronts 
and  major  disturbances  can  cause  large  variations  in  wind 
velocity  and  direction.  In  fact,  when  the  hotline  locator 
model  was  exercised  with  the  T+l  hr  Dean  and  Olmstead  winds, 
the  "wafer  ctrs"  hotline  did  move  northward.  Figure  25 
shows  the  hotline  predictions  assuming  the  T+l  hr  winds. 

To  correctly  represent  the  wind  structure,  a  model  would 
have  to  simulate  the  wind  profiles  all  along  the  path  of 
a  falling  particle.  Section  VII  describes  a  technique  to 
linearly  interpolate  temporal  wind  changes  in  model  layers 
beneath  falling  particles.  The  technique  was  used  to 
predict  the  segment  3  asymptote  with  a  linear  temporal 


Figure  25.  Hotlines,  Castle-Bravo,  T+lHr  Dean  and  Olmstead  Winds 


interpolation  of  shot-time  and  T+l  hr  wind  profiles  during 
particle  fall  and  translation. 

Sensitivities 

The  sensitivity  of  hotline  location  was  examined  for 
two  factors: 

1.  Pancake  cloud  altitude 

2.  Wind  profile. 

Pancake  cloud  height  variations  caused  consistent 
movement  of  the  hotline  to  directions  influenced  primarily 
by  the  mean  "effective"  wind  beneath  the  pancake.  Figure  26 
shows  the  hotline  from  a  pancake  cloud  located  at  30  km 
altitude.  The  agreement  with  the  actual  hotline  is  inter¬ 
esting,  but  it  is  unlikely  that  a  30  km  high  pancake  cloud 
could  have  generated  the  Bravo  fallout  hotline.  A  70  pm 
particle  takes  19.5  hours  to  fall  from  30  km  in  a  tropical 
atmosphere,  whereas  Figure  13  indicates  that  fallout  parti¬ 
cles  arrived  in  the  atolls  4  to  6  hours  after  detonation 
(Ref  14). 

Hotline  location  sensitivity  to  the  assumed  wind  pro¬ 
file  was  addressed  by  computing  hotlines  and  stack  traces 
with  different  winds.  Figures  27  and  28  show  that  hotline 
predictions  would  change  dramatically  if  the  particles  were 
translated  in  the  original  Three  Point  Average  winds  vs 
the  Dean  and  Olmstead  winds.  Three  point  average  winds  and 
components  are  plotted  in  Figure  29. 


Figure  26.  Hotlines,  Castle-Bravo,  30km  Pancake  Cloud 


Figure  27.  Hotlines,  Castle-Bravo,  Three  Point  Average  Winds 


RLTITUDE  (KM) 
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Results  of  the  Castle-Bravo  case  study  are  summarized 
as  follows: 

1.  The  hotline  locator  model  agreed  with  the  DELFIC 
prediction;  neither  predicted  the  hotline  curva¬ 
ture. 

2.  A  vertical  stack  technique  was  developed  to 
back  calculate  particle  starting  heights  in  the 
stabilized  cloud  to  approximate  the  actual 
hotline. 

3.  The  back  calculated  starting  heights  provided  a 
new  interpretation  of  the  stabilized  cloud  that 
included  a  stem  and  a  mixing  region  beneath  the 
cloud. 

4.  Hotline  location  is  extremely  sensitive  to 
the  assumed  winds  and  spatial/temporal 
variations . 


VII .  Linear  Temporal  Interpolation  for  Wind  Updates 

As  suggested  in  previous  sections,  the  use  of  a  single 
wind  profile  may  be  inadequate  for  predicting  fallout 
locations  at  distances  greater  than  the  synoptic  weather 
scale.  Wind  updating  can  be  achieved  with  a  pseudo- Lagrangian 
scheme  of  following  the  fallout  particles.  In  the  hotline 
locator  model,  a  particle  resides  in  a  layer  for  a  finite 
time.  During  that  time,  the  winds  in  layers  beneath  the 
particle  may  change.  If  the  changed  wind  profile  is  known 
at  a  later  time,  then  the  winds  in  any  layer  can  be  estimated 
by  linearly  interpolating  (in  time)  the  vector  component 
magnitudes . 

Specifically,  for  two  wind  profiles,  the  temporally 
interpolated  wind  speed  in  a  model  layer  is  computed  as 
follows : 

ULTI  =  + 

ULti  =  linear  temporal  interpolated  wind  speed 
UQ  =  wind  speed  at  time  zero 
Ut  =  wind  speed  at  time  t  >  0 

t  =  time  when  Ut  is  applicable  (if  t  ■  0 , 

there  is  only  one  wind  profile  and  the 
equation  does  not  apply) 
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T  -  elapsed  time  of  particle  fall  from  starting 
height  to  layer  where  interpolated  winds 
are  computed 

when  T  «  0  ,  UTTT  -  U 

Ll  I  O 

and  T  *  t  ,  ULTI  =  Ut 

The  technique  was  tested  with  the  T  =  0  and  T  +  1  hr 
Castle-Brave  (Dean  and  Olmstead)  winds  (Figure  30)  in  the 
hotline  locator  model.  In  Section  VI,  the  later  winds  were 
shown  . v  explain  Bravo  hotline  trend  beyond  280  km  from 
ground  zero.  This  exercise  was  an  attempt  to  check  the 
linear  temporal  interpolation  (LTI)  method  and  to  explain 
the  3rd  segment  of  the  AFSWP  hotline.  Stack  traces  are  shown 
in  Figure  31.  The  effect  of  linear  interpolation  to  the 
T  +  1  hr  wind  profile  is  that  particles  starting  from  high 
altitudes  land  to  the  north. 

The  hotlines  from  pancake  clouds  and  distributed  initial 
clouds  are  shown  in  Figure  32.  "LTI-P"  and  "LTI-wfrs”  refer 
to  LTI  winds  beneath  pancake  and  wafer-distributed  clouds, 
respectively.  The  "LTI-wfrs"  hotline  agrees  with  the  AFSWP 
hotline  beyond  280  km,  implying  that  particles  falling  from 
the  high  cloud  (not  stem)  were  more  affected  by  the  T  +  1  hr 
winds  than  by  the  T  +  0  winds.  However,  T  +  1  hr  winds  may 
not  be  recent  enough,  since  arrival  times  at  x  >  280  km  are 
much  longer  than  1  hour. 
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HDJOINT 

VERTICAL  STACKS 


Figure  31.  Ground  Traces  of  Single-Size-Particle  Stacks 

Dean  and  Olmstead  T+lKr  Winds 


Figure  32.  Hotlines,  Castle-Bravo,  Linear  Temporal  Interpola 

and  Fixed  Profile  Traces 


The  back  calculated  hotline  was  computed  by  combining 
the  linear-temporal  interpolation  technique  with  wafer  start¬ 
ing  heights  shown  in  Table  IV.  Figure  33  depicts  the  hotline 
predicted  with  the  the  locator  model  using  starting  heights 
estimated  from  stack  trace  intersections  and  linearly 
interpolated  T  +  0  and  T  +  1  hr  Dean  and  Olmstead  winds. 

The  calculation  properly  replicated  the  southward  extent  of 
the  hotline  observed  in  the  AFSWP  data.  The  wind  interpo¬ 
lation  produced  nearly  identical  hotlines  beyond  280  km 
from  ground  zero.  In  fact,  beyond  280  km,  particle  starting 
heights  were  predicted  with  the  correlations  derived  from 
DELFIC  data  (Section  III).  The  pancake  cloud  was  found  to 
best  represent  the  hotline  between  100  and  120  km  east  of 
the  burst  point.  Closer  than  100  km,  the  stack  trace  method 
indicated  that  the  stabilized  cloud  did  not  contain  gravity- 
sorted  activity  peaks.  In  fact,  analysis  of  the  stack  trace 
calculations  showed  that  some  large  particles  started  from 
higher  altitudes  than  smaller  particles.  The  low  altitudes 
in  which  this  occurred  correspond  to  the  stem  region  hypothe¬ 
sized  in  Section  VI. 

All  references  which  discussed  the  Castle-Bravo  test 
emphasized  the  large  uncertainties  in  the  fallout  trace  data; 
most  of  the  fallout  was  deposited  over  the  ocean  and  the 
measurement  netowrk  was  modest.  However,  the  southward 
curvature  was  verified  by  data  from  inhabited  atolls  (Rongerik 
and  Rongelap)  downwind  of  the  detonation.  At  close  range 
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Figure  33.  Hotlines,  Castle-Bravo,  Three-Segment  Fit 


(<  80  km),  it  is  interesting  to  note  that  the  stack  traces 
show  that  low  altitude  debris  would  first  move  to  the  south¬ 
west,  then  turn  to  the  east-northeast.  This  implies  that 
the  hotline  inferred  from  the  AFSWP  contours  (which  starts 
going  due  east  from  the  burst  point)  may  be  incorrectly 
located  and  should  actually  be  curving  southward. 
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VIII.  Solution  of  Dose  Rate  Integral  for  a  Curved  Hotline 

This  section  presents  the  second  step  to  modeling 
variable  winds  in  fallout  smearing  codes.  An  analytic 
expression  is  developed  for  the  activity  distribution  on 
the  ground  away  from  the  curved  hotline. 

The  fundamental  equation  for  ground  dose  rate  from  a 
smearing  code  was  given  as  Eq  (II-l). 


D(x,y) 


f(x,y,t')g(t')dt' 


(II-l) 


The  spatial  distribution  function  f(x,y,t')  was  given 
by  Eq  (II-2)  : 


f(x,y,t')  - 


SH  °x(t') 


EXP 


x-v  t' 


/TF  Oy,(t') 


EXP 


y-v  t' 


(11*2) 


f(x,y,t')  can  be  generalized  to  variable  winds  by 
replacing  the  terms  vxt'  and  v  t'  : 


v  t'  ->  /  v  dt' 


V' 


/  v  dt 
0  y 
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where  v  and  v  are  now  the  variable  wind  components 
x  y 

experienced  by  the  falling  particles  during  their  residence 
times  in  discrete  model  layers. 

Therefore,  for  a  curved  hotline  (v^  1  0),  the  spatial 
distribution  function  has  the  form: 


f(x»y»t)  = 


/7¥  a. 


EXP 


r  /  t 
fx-jv  dt 
0  x 


/Tir 


EXP 


t' 

/y-Jv  dt 
0  y 


(VIII-1) 


Reference  1  showed  that  g(t)  can  be  approximated  with 
the  first  two  terms  of  a  Taylor  series  expanded  about  t  , 
the  arrival  time  of  cloud  center  at  hotline  location  (X,Y)  : 

g(t-ta)  =  g(ta)  +  (||-)  (t-ta)  +  .  .  . 

La 

Then, 


D 


/  f(x-/vxdt~)f Cy-Jv  dt~)(jf0  (t 


ta)dt 


(VI II -2 ) 
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The  second  term  on  the  right-hand  side  of  (VIII-2)  is  zero 

CO  00 

if  the  limits  of  integration  are  changed  from  /  to  / 

0  -  00 

The  limit  change  is  valid  if  the  contribution  to  the 

00 

integral  from  the  /  portion  is  small.  The  contribution 

-00 

will  be  small  if  the  integral  is  evaluated  at  (x,y)  distant 
from  the  origin,  i.e.,  late  times  after  burst.  The  second 
integral  takes  the  form: 


00 

/  (gaussian) (gaussian) (constant ) (odd  function)dt  , 

-  CO 


which  is  the  integral  of  an  odd  function  about  the  origin, 
equal  to  zero. 

Only  the  first  integral  on  the  right-hand  side  of 
(VIII-2)  remains  to  contribute  to  the  dose  rate  calculation. 


D 


KYKgCta) 


00 


/ 

0 


1 

/2¥  ox 


EXP 


1 

/Z¥  oy 


EXP 


dt 


(V 1 1 1  -  3  ) 


Since  the  distribution  is  gaussian  centered  on  the 

t'  t' 

hotline,  the  displacements  /v  dt  and  /v  dt  are  coordinates 

0  x  0  y 

of  the  hotline  at  arrival  time,  t'  =  ta  .  These  coordinates 

4  •> 

*  > 
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were  determined  from  the  hotline  locator  model  (Section  IV) 
Therefore , 


a 

/  v  dt 
0  x 


X  =  V. 


r  vydt 

0 


Y  =  V  t 

y  a 


where  V  and  V  are  particle  residence  time  weighted 
x  y 

averages  of  wind  components. 


V 


x 


dX 

at 


Exponential  and  differential  terms  of  the  integrand 
become : 


EXP 


=  EXP 


Expanding  the  exponents  in  the  integrand: 


4  «• 


*  P 
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EXP  -li(jL2.-2xX>x2  + 

V  °x! 


V  v 

yl-2y(yl)x  +  (^)x2 
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i  . 
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=  EXP  -*f^T  ♦  HXP  X 
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x  +  x  X'  1  v: 

o  2  a  *  )  7  ( a  i  j-l 
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-I _ 

o  2a  2 

,  x  y 


EXP  X 


o  2x+a  2 
y  x 


V 
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a  2  o  2 
x  y 


X2/  y  ‘V 


,2+(rA)  a  2' 


a  2o  2 
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The  right  exponential  has  the  form:  EXP(-XB-X2A)  .  If 

oo 

the  contribution  from  /  is  negligible  (as  assumed  earlier, 

~  OO 

and  in  (1)),  the  integral  can  be  evaluated  as  one  of  the 


form: 


Jn  =  /  xn  EXP[-XB-X2A]dX 


with  n  =  0  . 
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The  solution  is  shown  in  Reference  16  as 


jo  -  VF  EXP[B2/4AD 
In  this  case: 


1 


Simplifying  the  exponential  term: 


[°y2x+ox2y(yZ)J4 


T 


2(^y2+ 
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The  integral  for  D  becomes: 


.  1  1^ 

/TiT  o  2V  2  +  o  2V  2 
y  x  x  y 


/Tif 


(a  2V  2+o  2V  2 
y  x  x  y 


-h 


(V 1 1 1  -  5 ) 


Combining  the  exponential  terrar  of  part  II  of.  the  integral: 
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Combining  the  constant  and  the  exponential  terms  for  D: 


D 


KV<V 

/Tn 


(o  2V  2+o  2V  z)'h 
y  x  x  y 


-is 

(xVy  -  yVx)2 

a  2V  2  +  a  2V  2 

> 

L  y  X  X  y 

\ 


(V 1 1 1  -  6  ) 


Equation  (VIII-6)  reduces  to  the  equation  described  in 
Reference  1  for  a  constant,  unidirectional  wind  when  V 


D  = 


KY  g(ta) 
/Tn  a  V 

y  * 


EXP 


cf*) 2 

y 


0: 


Reference  21  demonstrated  that  (VIII-6)  is  a  gaussian 
function  centered  at  the  hotline  coordinate  (x,y)  and 
perpendicular  to  a  straight  line  between  the  origin  (ground 
zero)  and  (x,y). 
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kyk  «<*„) 

/Tn 
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EXP 
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(xY  -  yX)2 

*2 

a  2X2  +  a  2Y2 
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_  y  x 

(VII I -7) 


Either  (VIII-6)  or  (VIII-7)  can  be  used  in  the  smearing 
codes  with  the  hotline  locator  model  to  compute  dose  rates 
at  the  ground  resulting  from  variable  winds. 
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IX.  Summary 

A  method  was  developed  to  enable  smearing  codes  to  predict 
fallout  footprints  in  variable  winds.  Two  steps  are  required: 
predict  the  curved  hotline  and  compute  the  off-axis  dose  rates. 

A  hotline  locator  model  was  developed  to  predict  curved 
hotlines. 

a.  The  hotline  locator  model  uses  real  wind  profiles 
and  the  McDonald-Davies  method  to  predict  particles' 
terminal  velocities  and  translations  through  atmospheric 
model  layers  to  the  ground.  Polynomial  fits  to  DELFIC 
data  are  used  to  compute* yield -dependent  slopes  and 
intercepts  of  a  linear  particle  size  vs  height  descrip¬ 
tion  of  the  stabilized  fallout  cloud.  The  assumption 
that  particle  heights  are  locations  of  constant-size 
activity  distribution  peaks  is  verified  by  results 
which  agree  with  DELFIC  output  (Washington  DC,  100  kt). 

b.  A  vertical  stack  technique  was  devised  to  perform 
adjoint  calculations  of  the  Castle-Bravo  hotline. 
Calculations  show  that  the  Bravo  hotline  can  be  explained 
by  the  presence  of  a  stem  and  a  mixing  region  beneath 
the  classical  fallout  cloud. 

c.  A  pseudo-Lagrangian  technique  was  developed  to 
update  atmospheric  winds  in  model  layers  beneath  the 
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falling  particles.  The  technique  was  tested  with 
T  =  0  and  T  +  1  hr  Bravo  winds  and  the  resulting  hot¬ 
line  agreed  with  the  late-time  AFSWP  hotline  asymptote. 


The  two-dimensional  dose-rate  integral  has  been  solved 
analytically,  enabling  smearing  codes  to  compute  off-axis 
levels  near  curved  hotlines. 


I 

I 


X.  Conclusions 


On  the  basis  of  the  research  and  calculations  performed 
for  this  study,  it  is  evident  that: 

a.  A  pancake  cloud  falling  through  an  altitude- varying 
wind  field  produces  a  slightly  curved  hotline.  Curvature 
is  minimal  and  does  not  improve  the  single,  effective 
wind  results  of  smearing  codes. 

b.  Realistic  hotline  curvature  can  be  obtained  by 
computing  the  motions  of  fallout  particle  arrays  which 
start  from  altitudes  that  are  the  geometric  centers  of 
the  DELFIC  population  of  wafers  in  stabilized  clouds. 
Starting  heights  can  be  predicted  from  polynomial  and 
least-squares  fits  to  DELFIC  data. 

c.  Particle  starting  heights  represent  altitudes  of 
peaks  in  constant- size  activity  distribution.  Therefore, 
the  ground  trace  of  the  activity  peaks  describes  the 
fallout  hotline. 

d.  Hotline  location  does  not  depend  on  particle  size 
distribution.  The  hotline  locator  model  prescribes  an 
array  of  trace  particles  only  to  define  the  hotline 
coordinates . 
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e.  Analysis  of  the  Castle-Bravo  fallout  data  and 
reconstruction  of  particle  trajectorires  inferred  that 
the  hotline  was  generated  from  (at  least)  three  distinct 
particle  size-altitude  functions:  stem,  mixing  region, 
and  stabilized  cloud. 

f.  Fallout  from  high  yield  shots  remains  aloft  so  long 
that  use  of  a  single  wind  profile  (much  less  a  single 
wind  vector)  may  introduce  significant  errors  into 
fallout  predictions.  A  linear  temporal  interpolation 
technique  has  been  developed  to  update  winds  beneath 
particles,  as  the  particles  fall  and  translate. 

g.  The  two-dimensional  dose  rate  integral  was  solved 
analytically  to  produce  a  form  which  reduces  to  the 
one-dimensional  form  when  a  constant  wind  vector  is 
assumed. 

h.  The  two-step  method  will  enable  smearing  codes  to 
predict  realistically  curved  hotlines  with  real  wind 
data  and  to  compute  dose  rates  on  and  near  the  hotline. 
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Appendix  A 


Atmospheric  Models 

Atmospheric  state  parameters  are  used  to  compute  terminal 
velocities  of  fallout  particles  at  appropriate  heights  in  the 
atmosphere.  In  this  study,  the  computational  equations  from 
Reference  9  were  used  to  generate  altitude-dependent  density 
and  kinematic  viscosity  values  needed  to  calculate  particle 
fall  speeds.  Following  is  a  description  of -the  process  used 
to  compute  atmospheric  terms  in  the  U.S.  Standard  and  a  tropi¬ 
cal  atmosphere. 

U.S.  Standard  Atmosphere 

Altitude-dependence  of  the  atmospheric  state  variables 
is  motivated  by  the  assumed  temperature  profile.  The  standard 
atmosphere  uses  the  temperature  profile  shown  in  Figure  34; 
temperature  decreases  linearly  from  the  earth's  surface  to  a 
geopotential  height  of  11  kilometers,  at  which  point  the  tem¬ 
perature  profile  is  isothermal  up  to  20  kilometers.  The 
computational  sequence  is  as  follows: 


1.  Geopotential  Height 


HGEO 


r*Re*ZMSL 

(Re+fcMSL) 


(A- 1 ) 
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HGEO 


r 

Re 


geopotential  height  (m) 

acceleration  of  gravity  in  geopotential  meters 
acceleration  of  gravity  in  meters 


*0 

effective  earth  radius  a  6356766  meters 


ZMSL  =  height  of  calculation  (meters  above  mean  sea 
level) 


2 .  Temperature 

TZ  =  T  .  +  L  .  (HGEO-HB)  (A- 2) 

m ,  d  m ,  b 

TZ  =  temperature  at  ZMSL  (0°K) 

T  =  temperature  at  base  of  profile  segment 
m,b 

=  288.15  (°K)  for  HGEO  5  11  km 

=  temperature  at  11  km  for  11  km  <  HGEO  -  20  km 

0  K 

L  =  lapse  rate  in  profile  segment  (me~ter) 
m,b 

-  -.0065  —  for  HGEO  <  11  km 

m 

=  0  for  11  <  HGEO  5  20  km 

HB  «  geopotential  height  of  temperature 
profile  segment  base 

-  0  meters  for  HGEO  -  11  km 

■  11000  meters  for  11  km  <  HGEO  *  20  km 
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3.  Pressure 


00  If  Lm_b  f  0, 


PR  =  PB 


•  [T^r^tWTnrr] 


[8o'Mo  1 


(A-  3) 


PR  *  pressure  at  AMSL  (Newtons/m2  5  Pascals) 


pressure  at  base  of  temperature 
profile  segment  (Pascals) 


9  80665  (ge°P°tential  meters^ 
second2 

molecular  weight  of  air  at  sea  level 

-  28-9644  raft* 


constant  =  8.31432x10  (VN  l  ?») 


(b)  If  Lm>b  »  0, 


PR  =  PB 


*  EXP 


g0'M0(GHE°-HB) 


4.  Density 


DENS 


PR  •  M. 


DENS  ■  Mass  density  (—**■) 

m3 


(A-4) 


(A- 5) 
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5 .  Kinematic  Viscosity 


VISK 


VISK 

B 

S 


(B  (TZ) 1  *  5  /  ncvie 

(Tz+sy _  j  DENS 

2 

Kinematic  viscosity  (“— ) 

5  vv 

constant  =  l,458xl0'6  (kg/ (sec ‘m*0^) ) 
Sutherland  constant  =  110°  K 


(A-6) 


Tropical  Atmosphere 

Most  of  the  high-yield  atmospheric  test  data  was  obtained 
from  nuclear  weapon  tests  in  the  Bikini  Atoll.  Bikini  is 
located  at  11.35°  north  latitude  and  165.20°  east  longitude 
(Ref  10).  To  compute  particle  motions  in  realistic  atmospheric 
parameters,  a  tropical  atmosphere  was  computed  using  the 
equations  from  Reference  9.  However,  a  tropical  temperature 
lapse  rate  was  used  to  motivate  the  altitude-dependence  of  the 
state  variables.  Specifically,  TZ  is  obtained  from  the  seg¬ 
mented  temperature  profile  given  in  Reference  11  and  shown 
in  Table  A- I.  Figure  34  shows  the  tropical  lapse  rate  compared 
to  the  U.S.  Standard  profile.  The  lapse  rate  reported  for  the 
Castle-Bravo  Pacific  test  is  also  plotted  on  Figure  34  to  show 
the  close  agreement  between  computed  and  actual  values. 

Figures  35,  36,  and  37  plot  pressure,  density  and  kinematic 
viscosity  computed  for  both  atmospheres.  Figure  35  includes  the 
pressure  profile  reported  for  the  Castle-Brave  weapon  test. 
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Table  A-I 


Temperature  vs  Altitude:  Tropical  Atmosphere 


Altitude  (km) 

Temperature 

(°C) 

0 

+  26.5 

2.256 

+  13.0 

2.507 

+  13.8 

16.582 

O 

o 

00 

1 

22.128 

1 

VI 

oo 

o 

47.463 

-  3.0 

A  vertical  temperature  profile  is  generated  by 
linear  interpolation. 
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DENSITY  (KG/M3) 


Figure  36.  Density  vs  Altitude 
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Figure  37.  Kinematic  Viscosity  vs  Altitude 
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Appendix  B 


Hotline  Locator  Model  Program 

The  hotline  locator  model  was  written  in  FORTRAN  V 
and  executed  from  the  AFIT  computer  facility  using  the  ASD 
CDC  6600  computer.  Ambient  wind  data  is  read  into  the 
code  by  defining  the  reported  wind  speed  and  direction  at 
integer  levels  of  the  model  (0  to  30  km  in  100  meter 
increments).  See  lines  140  through  196.  This  listing 
contains  the  Dean  and  Olmstead  wind  data  developed  for 
the  Castle-Bravo  nuclear  test  (Ref  4).  This  appendix  con¬ 
tains  : 

I.  The  computational  algorithm 

II.  A  list  of  code  variables  with  definitions 

III.  A  listing  of  the  main  program. 
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I.  The  Computational  Algorithm 

The  hotline  locator  model  uses  real  wind  data,  a  twenty 
particle  size  array  and  the  McDonald-Davies  method  to  calculate 
ground  traces  of  particles  which  start  from  specific  altitudes 
in  a  stabilized  nuclear  cloud.  As  contained  in  Part  III,  the 
model  computes  according  to  the  following  algorithm: 

1.  Read  Inputs 

a.  Particle  density  (g/cm3) 

b.  Twenty  particle  diameters  (cm) 

c.  Actual  hotline  coordinates  (x,y  in  kilometers 
east  and  north  of  ground  zero) 

d.  Pancake  cloud  levels  (three  model  layer 
numbers  5  300;  high/low/increment;  hotline 
from  a  single  pancake  cloud  will  be  computed 
when  low  +  increment  >  high) 

e.  Twenty  particle  starting  levels  for  adjoint 
calculation  (model  layer  numbers  specify  the 
altitudes  from  which  each  particle  starts 

to  fall) 

f.  Weapon  yield  (kilotons) 


2.  Compute  atmospheric  density  and  kinematic  viscosity 
from  equations  described  in  Section  III. 


Read  actual  wind  velocities  (m/sec)  and  directions 
(degrees  clockwise  from  north)  at  reported  atmospheric 
heights.  Velocities  and  directions  are  entered  at 
appropriate  model  levels,  where  levels  are  separated 
by  100  m  and  numbered  from  1  to  301. 

Compute  x  and  y  components  of  the  wind  vector. 

Since  standard  meteorological  practice  is  to  report 
wind  direction  in  "<|>  degrees  clockwise  from  north- 
upwind,"  the  vector  components  of  wind  are: 

Vx  =  V  (-sin  <{>) 

Vy  »  V  (-cos  <}>) 

Linearly  interpolate  the  component  magnitudes  between 
reported  levels  and  assign  two  components  to  each 
model  level.  Compute  the  interpolated  wind  angles 
at  each  level  to  compare  with  actual  profiles. 

Use  Davies -McDonald  to  compute  terminal  velocities 
of  each  particle  size  at  each  model  level. 

Compute  residence  time  of  each  particle  in  each 
layer. 

Compute  the  two-dimensional  translation  of  each 
particle  in  each  layer. 


9.  Compute  residence- time  weighted  wind  vectors  for 
comparison  with  the  method  described  in  References 
2  and  8.  The  effective  wind  vector,  for  uni¬ 
directional  smearing  code  calculations,  is  obtained 
by  "weighting"  the  altitude-distributed  winds  by 
the  time  that  particles  reside  in  the  corresponding 
altitudes.  Since  particles  fall  faster  at  higher 
altitudes,  high  altitude  winds  contribute  less  to 
the  effective  wind  vector  than  low  altitude  winds. 
The  effective  wind  is  a  mean  value  defined  as: 


n 


l  t.W 


i  =  l 


l  l 


T  . 
X 


W  =  mean  effective  wind,  speed  (m/sec) 

=  residence  time  of  a  particle  in  i**1 
layer  (second) 

*  wind  speed  in  i**1  layer 


10.  Initialize  Plot  #1  by  calling  plot  library 
subroutines . 


11.  Compute  and  plot  hotline  coordinates  for  particles 
dropped  from  pancake  cloud  level  through  a  real  wind 
distribution.  Hotline  coordinates  are  obtained  by 
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summing  the  x  and  y  translations  in  all  model  layers 
beneath  the  specified  pancake  cloud  level. 

12.  Compute  slope  and  intercept  of  particle-height  lines 
from  polynomial  fits  with  specified  yield. 

13.  Use  slope  and  intercept  to  compute  20  starting 
altitudes  of  20  particle  sizes  in  stabilized  clouds. 

14.  Compute  and  plot  hotline  coordinates  for  particles 
which  fall  from  heights  computed  in  Step  13. 

15.  Plot  the  real  hotline  coordinates  provided  as  input. 

16.  Initialize  Plot  #2. 

17.  Plot  the  real  hotline  coordinates  provided  as  input. 

18.  Compute  ground  traces  for  each  of  20  particles  dropped 
through  each  of  the  300  model  layers. 

19.  Plot  the  ground  traces  of  vertical  stacks  of  particles 
which  extend  from  the  surface  to  adjoint  heights 
(input  #l.e).  Vertical  stacks  produce  ground  contours 
which  help  to  define  the  altitudes  from  which  particles 
would  have  to  fall  in  order  to  land  on  the  real 
hotline . 

The  model  is  designed  as  a  diagnostic  program  to  compute  and 
plot  hotline  coordinates  and  stack  traces  after  printing 
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intermediate  results,  such  as  wind  profiles,  particle 
velocities  and  atmospheric  state  variables. 
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II.  Code  Variables  and  Definitions 


1. 

A 

"** 

The  number  of  atmosphere  layers  in  a  vertical 
stack  (back  calculations) 

2. 

ADXT 

ss 

Two  dimensional  array  of  ground  coordinates 
(x-direction)  for  each  particle  size  and  model 
level  at  starting  height 

3. 

ADYT 

Two  dimensional  array  of  ground  coordinates 
(y-direction)  for  each  particle  size  and  model 
level  at  starting  height 

4. 

AFX 

= 

x-coordinate  array  of  the  actual  fallout  hotline 

5. 

AFY 

= 

y-coordinate  array  of  the  actual  fallout  hotline 

6. 

AP 

= 

Array  of  model  levels  that  specify  adjoint 
stack  heights 

7. 

BETA 

= 

Constant  used  to  compute  kinematic  viscosity 
of  air 

8. 

DENS 

= 

Array  of  atmospheric  density  values  at  each 
of  301  model  levels 

9. 

DVX 

= 

Incremental  velocity  in  x-direction  computed 
between  published  wind  levels  and  used  for 
linear  interpolation  fallout  wind  profiles 

10. 

DVY 

a 

Incremental  velocity  in  y-direction  computed 
between  published  wind  levels  and  used  for 
linear  interpolation  fallout  wind  profiles 

11. 

DX 

= 

Two  dimensional  array  of  x-distances  traveled 
by  each  particle  in  each  model  layer 

12. 

DXT 

a 

Array  of  ground  coordinates  (x-direction) 
obtained  by  summing  DX  from  the  ground 
upward 

13. 

DY 

a 

Two  dimensional  array  of  y-distances  traveled 
by  each  particle  in  each  model  layer 

14. 

DYT 

= 

Array  of  ground  coordinates  (y-direction) 
obtained  by  summing  DY  from  the  ground  upward 

15. 

GAMMA 

a 

Ratio  of  gravitational  constants:  meters  to 
geopotential  meters 
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16, 

GOP 

* 

Graviational  constant  in  geopotential  meters 
per  second  squared 

17. 

HB 

“ 

Bottom  height  of  an  atmospheric  temperature 
profile  segment 

18. 

HGEO 

■ 

Geopotential  height  in  atmosphere 

19. 

I 

* 

Dummy  index  for  manipulating  subscripted 
variables 

20. 

IL 

a 

Natural  logarithm  of  intercept  obtained  from 
polynomial  least  squares  correlation  with 
natural  logarithm  of  weapon  yield 

21. 

INTCPT 

* 

Height  of  zero  diameter  particle  in  a 
stabilized  cloud:  EXP(IL) 

22. 

IPAK 

3 

Array  of  characters  used  by  DISSPLA  to  store 
information  for  plot  legends 

23. 

J 

= 

Dummy  index  for  manipulating  subscripted 
variables 

24. 

L 

S 

Dummy  index  for  manipulating  subscripted 
variables 

25. 

LMB 

s 

Slope  of  an  atmospheric  temperature  profile 
in  °K/meter 

26. 

MO 

= 

Mean  molecular  weight  of  air 

27. 

N 

s 

Dummy  index  for  manipulating  subscripted 
variables 

28. 

PB 

z 

Atmospheric  pressure  (Newtons/m2)  at  the  base 
of  each  linear  temperature  profile  segment 

29. 

PCB 

s 

Input  height  (model  level)  of  a  pancake  cloud: 
the  lowest  level  desired  for  calculation  of 
hotline  coordinates 

30. 

PCS 

s 

Input  incremental  height  (model  levels)  between 
pancake  clouds,  if  multiple  calculations  are 
computed 

31. 

PCT 

s 

Input  height  (model  level)  of  a  pancake  cloud: 
the  highest  level  desired  for  calculation  of 
hotline  coordinates.  Note:  If  PCS  >  (PCT-PCB) 
the  code  will  compute  the  hotline  coordinates 
from  a  single  pancake  cloud  (PCB). 
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32.  PDEN 


Fallout  particle  density  (grams/cm3) 

Array  of  20  fallout  particle  diameters  (cm) 


33. 

PDIA 

s 

Array  of  20  fallout  particle  diameters  (cm) 

34. 

PDXT 

St 

Array  of  hotline  x-coordinates  used  for  plot 
subroutine 

35. 

PDYT 

= 

Array  of  hotline  y-coordinates  used  for  plot 
subroutine 

36. 

PHI 

= 

Array  of  wind  directions  in  degrees  clockwise 
from  north:  upwind 

37. 

PR 

= 

Array  atmospheric  pressure  values  (Nt/m2) 
at  each  model  level 

38. 

Q 

= 

Product  of  particle  drag  coefficient  and  the 
square  of  Reynold’s  number  used  for  terminal 
velocity  calculation  with  Davies'  method 

39. 

R 

= 

Array  of  Reynold’s  numbers  computed  with 

Davies'  correlations  with  Q 

40. 

RERTH 

= 

Earth  radius  (6356766  meters) 

41. 

RL 

— 

Array  of  (log^R)  computed  with  Davies' 
correlations  with  log10Q 

42. 

RSTAR 

s 

Universal  gas  constant  (8314.32(N+  M)/ (KMole*°K) 

43. 

RTD 

= 

Array  of  radial  distances  from  burst  ground 
zero  to  ground  impact  coordinates  of  each 
particle  size 

44. 

RX 

= 

X-coOrdinate  array  of  the  actual  fallout  hotline 

45. 

RY 

= 

Y-coordinate  array  of  the  actual  fallout  hotline 

46. 

S 

= 

Sutherland's  constant  (110°  K)  used  to  compute 
kinematic  viscosity  of  air 

47. 

SL 

a 

Natural  logarithm  of  slope  obtained  from 
polynomial  least  squares  correlation  with 
natural  logarithm  of  weapon  yield 

48. 

SLOPE 

3 

Height  variation  of  particle  sizes  in  a 
stabilized  cloud:  EXP(SL) 

49. 

SVX 

= 

Array  of  x-components  of  wind  summed  to 
compute  a  mean  wind  (used  for  reference  only) 
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50. 

SVY 

SB 

Array  of  y-components  of  wind  summed  to 
compute  a  mean  wind  (used  for  reference 
only) 

51. 

TAU 

S 

Two-dimensional  array  of  residence  times  of 
particles  between  each  pair  of  model 
levels  (sec) 

52. 

TB 

= 

Temperature  at  the  base  of  atmospheric 
temperature  profile  segment  (°K) 

53. 

TZ 

= 

Array  of  atmospheric  temperature  values  at 
each  model  level 

54. 

VISK 

s 

Array  of  kinematic  viscosity  values 
(sec2/meter)  at  each  model  level 

55. 

VT 

Two-dimensional  array  of  terminal  velocities 
(cm/sec)  of  each  particle  size  at  each  model 
level 

56. 

VTA 

= 

Array  of  average  terminal  velocity  in  a 
model  layer  (cm/sec) 

57. 

VX 

= 

Array  of  wind  velocity  x-components  (m/sec) 

58. 

VXEFF 

Two-dimensional  array  of  effective  Vx  values 
obtained  by  averaging  the  residence-time- 
weighted  VX  values  from  each  model  layer 
(used  for  reference  only) 

59. 

VXM 

= 

Array  of  mean  wind  x-components  beneath 
each  model  level  (used  for  reference  only) 

60. 

VY 

= 

Array  of  wind  velocity  y-components  (m/sec) 

61. 

VYEFF 

3= 

Two-dimensional  array  of  effective  VY  values 
obtained  by  averaging  the  residence-time- 
weighted  VX  values  for  each  model  layer 

62. 

VYM 

= 

Array  of  mean  wind  x-components  beneath 
each  model  level  (used  for  reference  only) 

63. 

VZ 

*“ 

Array  of  terminal  velocities  of  each  particle 
size  obtained  from  definition  of  Reynold’s 
number  in  Davies'  method  (cm/sec) 

64. 

WIND 

* 

Array  of  wind  speeds:  magnitudes  of  wind 
vectors  at  each  model  level 

65. 

WT 

a 

Array  of  particle  weights  (g  cm/sec2) 

107 


66. 

WVEFF 

Two-dimensional  array  of  effective  wind 
speeds  obtained  by  averaging  the  residence¬ 
time-weighted  VX  values  for  each  model  layer 

67. 

WVM 

= 

Array  of  mean  wind  speeds  beneath  each  model 
level  (used  for  reference  only) 

68. 

X 

= 

Array  of  hotline  x-coordinates  used  for 
plotting  adjoint  stacks  (km) 

69. 

Y 

= 

Array  of  hotline  y-coordinates  used  for 
plotting  adjoint  stacks  (km) 

70. 

YKT 

= 

Weapon  yield  in  kilotons 

71. 

YL 

= 

Natural  logarithm  of  YKT 

72. 

Z 

Height  of  particle  size  in  a  stabilized 
cloud  generated  with  linear  fits  to  mean 
wafer  center  heights  and  polynomial  fits 
for  slope  and  intercept  (m) 

73. 

ZMSL 

= 

Array  of  model  level  heights  above  sea 
level  used  to  generate  atmospheric  parameters 
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